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We propose a new approach to quantum phase transitions in terms of the density-functional fi- 
delity, which measures the similarity between density distributions of two ground states in parameter 
space. The key feature of the approach, as we will show, is that the density-functional fidelity can 
be measured easily in experiments. Both the validity and versatility of the approach are checked by 
the Lipkin-Meshkov-Glick model and the one-dimensional Hubbard model. 
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Motivation — Rece ntly , considerable attentions llU2 J3L 
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have been paid to the role of fidelity [2fJ, [2l|, a con- 
cept emerging from quantum information theory [22| . in 
quantum phase transitions [23| . Since the fidelity mea- 
sures similarity of two quantum states without consid- 
ering any order information in the states, its approach 
to quantum phase transitions has shed new lights on 
our understanding of the quantum criticality. For exam- 
ple, the fidelity and its lea ding term show singular and 
scaling behavior [l3l EH EH, Il6j | in the topological phase 
transition 24( occurring in the ground states of both the 
Kitaev honeycomb model (2f| and the Kitaev toric model 
[26| . While the topoloigical phase transitions cannot 
be described by the Landau- Ginzburg- Wilson sponta- 
neous symmetry-breaking theory. Mathematically, the 
fidelity is defined as an overlap between two ground states 
[20L l2l| separated by a certain distance in the parame- 
ter space. To calculate the fidelity, one needs to know, 
in principle, all information of the ground states, in- 
cluding their each component in the Hilbert space of 
the Hamiltonian. This requirement is not an insuper- 
able barrier for theoretical studies because one can use 
some smart methods, for instance, the density-matrix- 
renormalization-group technique H3, HH, to reduce 
the dimension of the Hilbert space. However, it might be- 
come a serious problem for experimental studies because 
to obtain the full ground-state function of a quantum 
many-body system is extremely difficult in experiments. 
This fact, to a certainty, obstructs experimentalists from 
measuring the ground-state fidelity of realistic quantum 
many-body systems. 

Therefore, to find an appropriate approach that not 
only can describe the change in the ground-state struc- 
ture in perspective of information theory, but also is 
based on experimentally measurable quantities, becomes 
a nontrivial problem. In this paper, we propose a new 
approach to quantum phase transitions in terms of the 
density-functional [1^, [3(| fidelity (DFF) and its sus- 
ceptibility (DFFS). According to the Hohenberg-Kohn 
theorems [29|, the ground-state properties of a quantum 



many-body system are uniquely determined by the den- 
sity distribution n x that minimizes the functional for the 
ground-state energy Bojiix]. Therefore, the distribution 
n x captures the most relevant information of the ground 
state. Any change in the structure of the wavefunction 
can be found by calculating the similarity between two 
density distributions, i.e. the DFF, directly. Since the 
density distribution is usually measurable in experiments, 
our approach, therefore, provides a practicable strategy 
to study quantum critical phenomena both theoretically 
and experimentally. To check both the validity and ver- 
satility of the approach, we take the Lipkin-Meshkov- 
Glick (LMG) model [3~D ] and the one-dimensional Hub- 
bard model [32[ as examples, and study the DFF of real- 
space and momentum-space density distributions in the 
two models, respectively. We show that the real-space 
DFF is able to witness the second-order quantum phase 
transition [III, [33|, [HI occurring in the ground state of the 
LMG model in the anisotropic case; and the momentum- 
space DFF witness the Mott- insulator transition [3oT ]. 
which is of Beresinskii-Kosterlitz-Thouless (BKT) type 
[3^ . H3] , in the ground state of the Hubbard model. 

Formalism — To begin with, we consider a general 
Hamiltonian of quantum many-body systems 
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where Hj is the interaction term and A denotes strength, 
and n x is the local (pseudo) potential associated with 
density distribution {n x }. The index x can be discrete 
or continuous depending on the system under study. 
Though in the local-density-approximation (LDA) cal- 
culation, n x usually refers to the density of electrons in 
real space, it can also be generalized to population in con- 
figuration space of a reduced-density matrix or the den- 
sity of state in energy(momentum) space. The essence of 
the density-functional theory is that the density distribu- 
tion has already capture the most relevant information 
of the ground state, any change in the environment con- 
ditions will leads another unique density distribution. In 
quantum information perspective, the density-functional 
theory can be used to calculate the ground-state entan- 
glement in quantum many-body systems (38j. 

Mathematically, the single-particle state can be ob- 
tained by tracing out all other particles' degree of free- 
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dom, i.e. pixux'i) =tr|* (A))(* (A)|, where |* (A)> 
is the ground state of the Hamiltonian |T]). Then the 
density distribution n is simply the diagonal part of the 
reduced-density matrix p(x\,x' 1 ) 



\x)(x\, 



p(x,x) 



(2) 



where n has already been normalized. According to 
the Hellmann-Feymann theorem, the density distribution 
can also be obtained as 

n x = (#o(A)|nx|*o(A)) = (M^WH /dfi x )\M^))- 

(3) 

Therefore, the change in the ground-state wavefunction 
can be reflected from the fidelity between the two density 
distributions. For two ground states at A and A', the DFF 
has the form, 



F(X, A') = tTy/n(X)n(X'). 



(4) 



If we fix the distance 5X = A — A', the DFF is expected to 
show a drop around the critical point because two ground 
states in different quantum phases has the maximum dis- 
tance. In experiments, the density distribution n x (X) is 
measurable. Then the DFF between the density distri- 
butions under different environmental condition can be 
obtained. 

Meanwhile, it was found that the fidelity susceptibility 
which denotes the leading term of the fidelity, plays 
a central role in the fidelity approach to quantum phase 
transitions. Expanding the DFF to the leading order, we 
can find that 



F(X,X + SX) = 1 



(sxy 



-XF- 



where the DFFS xf takes the form 



XF 



E- 



dn x 
~~dX 



(5) 



(6) 



Therefore, if we regard dn x /dX as an independent func- 
tion besides the density distribution n x , the DFF is a 
functional of n x and dn x /dX, both of which, in prin- 
ciple, maximize the DFFS at the critical point. The 
typical case might be that the denisty n x in a certain 
region x m - m < x < x max vanishes (for l/n x ) rapidly (for 
dn x /dX). These conclude the main formulism of the DFF 
approach to quantum phase transitions. 

The Lipkin-Meshkov-Glick (LMG) model — To con- 
vince oneself of the validity of the approach, let us first 
study the DFF in a simple quantum phase transition oc- 
curring in the LMG model, which can provide us a clear 
paradigm of the approach. The Hamiltonian of the LMG 
model can be written as 
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FIG. 1: (color online) The density-distribution properties of 
the LMG model are analyzed for a sample of S = 1024 and 
7 = 0.5. (a) The density distribution n z as a function of s z 
for various h = 0.8, 0.825, 1.1 (from left to right). The 
rightmost peak is for h = 0.975 (b) The DFF as a function of 
h for Sh = 0.005. (c) The DFFS as a function of h. 



where S a (a = x,y,z) is the spin-S* operator and S± = 
S x ± iS y . The prefactor 1/5 is necessary to ensure a 
finite E/S in the large 5 limit. It is understood that 
the total spin is the conserved quantities, i.e.,\H, S 2 ] = 
0. The ground state of the Hamiltonian locates in the 
subspace of maximum S and can be expressed in the basis 
of S z , \^o(h)) = ^2f(s z )\s z ). Therefore, the ground- 
state energy is a density functional i£o[n«], where n z = 
(*oO)N*oW) with s z = \s z ){s z \. 

The ground state of the LMG model in the anisotropic 
case (7 7^ I) consists of two phases, i.e the polarized 
phase for h > 1, and the symmetry-breaking phase for 
< h < 1. A simple quantum phase transition occurs at 
the critical point h c = 1. For an illustrative purpose, we 
show the density distributions n z as a function of s z for a 
sample of S — 2048 under various hs in Fig. [2(a). When 
h < 1, the density distribution n z always shows a peak 
which is dragged to the boundary of maximum s z little 
by little as h increases. The DFF can be sketched out 
from the overlap between two neighboring distributions, 
which becomes smaller and smaller, and finally reaches a 
minimum around h c = 1. When h > 1, the density peak 
diverges. In this case, a small change of h can not change 
the density distribution significantly, then the DFF is 
close to 1 again. Fig. [1] (b) represents the DFF as a 
function of h for 8h = 0.005. Clearly the sudden drop of 
the DFF corresponds the second-order quantum phase 
transition. The sharp peak in the DFF susceptibility 
[Fig. Q] (c)] implies that the density distribution evolves 
dramatically around the critical point. 

For the LMG model, n z = ip* (s z )ip(s z ), so the DFF 
is mathematically the same as the previously studied 
ground-state fidelity [ll[, hence satisfies the same scal- 
ing and critical behaviors as the ground-state fidelity. 



However, we would like to emphsis here that, the mo- 
tivation of the DFF is quite different. The mathematical 
coincidence of two fidelities of the LMG model is due to 
that the model in the subspace of maximum S becomes 
a single-particle problem. In this case, the diagonal ele- 
ment of the single-particle reduced-density matrix is just 
the absolute value of the wavefunction. 

The one- dimensional Hubbard model — The density 
distribution in the density-functional theory can be gen- 
eralized to any density distribution, for instance, the den- 
sity of state. To see the versatility of the approach, here 
we take the one-dimensional Hubbard model [13, Hl| to 
show how the change in the density of state in momen- 
tum space around the critical point can be reflected from 
the corresponding DFF and DFFS. The Hamiltonian of 
the one-dimensional Hubbard model reads 



l,aCj>) + ^E%T«J,i ( 8 ) 



where a and c J;CT , a =|, J, are creation and annihilation 
operators for fermionic atoms with spin a at site j respec- 
tively, fiu = c^Co-, and U denotes the strength of on-site 
interaction. The Hubbard model can be solved exactly 
via the Bethe-ansatz method [35| . For the case of U > 
and N < L, the energy spectra of the system under the 
periodic boundary conditions are determined by a set of 
charge and spin rapidities {kj, A Q }, which satisfy the fol- 
lowing transcendental Bethe-ansatz equations (BAEs) 



M 



2nlj = kjL — #i(A a — sinfej), 



(9) 



N M 

2irJ a = J2°i ( A - - sin k i) ~ 9 ^ Xa ~ Afc )' ( 10 ) 

3=1 6=1 

where 8 n (k) = 2 tan _1 (4fc/n[7), M is the number of down 
spins and {Ij, J a } play the role of quantum number. The 
ground-state solution only consists of real ks and real As. 
It is a singlet state, given by a quantum number configu- 
ration (successive integers or half-odd-integers) symmet- 
rically arranged around zero. Once the rapidities {kj, A a } 
are obtained from Eqs. (|9I10[) . the energy of the state can 
be calculated as E — —2 J2j cos kj. 

For the half-filled case, i.e N = L, a Mott-insulator 
transition occurs at the critical point [7 = 0. If U > 0, 
a charge gap opens and the system becomes an insula- 
tor. The ground-state energy the Hubbard model around 
U = is infinitely differentiable, the phase transition be- 
longs to the BKT type. Recent entanglement approach 
to the Hubbard model [H, [4f| shows that single-site en- 
tanglement shows a maximum, due to the various corre- 
lations have long-range behaviors, at the critical point. 
However, whether the ground-state fidelity can witness 
the Mott-insulator transition is still controversial [f| . 
To apply the DFF approach to the phase transition, we 
use the concept of the density of state in the quasi- 
momentum space that is defined as p[(fcj_|_i + kj)/2] = 
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FIG. 2: (color online) The density of state and DFF of the 
Hubbard model are analyzed for a finite sample of L — N = 
210. (a) The density of state for various on-site interaction 
U = 20, 19.6, . . . , 0.04 with step 0.04. Two ends of each line 
denotes the momentum cut-off in finite systems. (b)The DFF 
as a function of U. (c) The DFFS a function of U. 



l/[L(kj + i — kj)]. Then the ground-state energy can be 
calculated as E = —(N/w) dkp(k) cos k. Clearly, the 
density of state p(k) uniquely decides the ground-state 
energy. 

To see the evolution in the density of state, in Fig. 
[2] (a), we show the density of state for a system of 
L = N = 210 and a set of uniformly distributed U with 
space SU = 0.04 . From the figure, we can see that if U 
is large, the densities of state are very close, this observa- 
tion means that the DFF between two closing densities 
is almost 1 [See Fig. [5] (b)]. The underlying physical 
picture is that, in the large U limit, the system becomes 
a full-filled spinless fermion model, the ks in Eq. are 
almost uniformly distributed, then the density of state is 
very flat in the momentum space. On the other hand, 
as U becomes smaller and smaller, the second term of 
Eq. ([9]) makes any k more close to the zero point hence 
suppresses the density of state in the region of \k\ > n/2, 
while increase p(k) in |fc| < 7r/2. The physics is straight- 
forward. The Pauli's exclusion principle does not for- 
bidden two electrons with different spin polarizations to 
occupy the same state in k space. To have a low en- 
ergy, electrons want to stay as close as possible to the 
k = point. Then the density of state become higher 
in |fc| < 7r/2 and smaller in |fc| > it/ 2 as U deceases. 
Around the critical point, the change in the density of 
state becomes dramatic. Then the DFF deviates from 
1 little by little, and the fidelity susceptibility becomes 
larger and larger. Though the Bethe-ansatz equations 
are extremely difficult to be solved around U = 0, we 
can still find the DFFS tends to a maximum as U tends 
to zero. 

We can also check the above picture analytically in 
some special cases. In the thermodynamic limit and half- 
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filled case, the density of state is [35| 

1 cosfc f°° dpJ {p)cos(psmk) 

p{k) = ^ + — J Q — TT^^ — ' 

where Jo is zeroth order Bessel function. If U = oo, 
the second term of Eq. (fTTj) vanishes, so p(k) = l/2n 
for k £ [— 7r, 7r], If J7 = 0, since J °° dpJo(p) cos(psin fc) = 
cosfcj -1 , so p(fc) — 1/tt for fc G [— 7r/2, tt/2] and ,o(fc) = 
for |fc| > 7r/2. This fact tells us that the density of state 
undergoes a dramatic change around the critical point. 
Therefore, as implied from the Eq. ([6]), the DFF should 
reaches a maximum at the critical point U = 0. 

In summary, we have proposed a new approach to 
quantum phase transitions based on the DFF and DFFS. 
The former measures the similarity between two density 
distributions and the latter describes the changing rate 
of the density distribution in the parameter space. We 
show that the divergence of the DFFS is typically related 
to a quickly vanishing density distribution, which actu- 



ally denotes a quantum phase transition. To check the 
validity of the approach, we studied the DFF in the LMG 
model and show that the divergence of the DFF at the 
critical point. Also to see the versatility of the approach, 
we studied the DFF for the density of state of the one- 
dimensional Hubbard model. As the on-site U tends to 
zero, the density of state in the regions [— 7r, — tt/2] and 
[7r/2,7r] quickly vanishes. This phenomena leads a maxi- 
mized (at least) DFFS around the critical point U = 0. 

Though we restricted out studies in two well-studied 
strongly correlated systems, the approach we proposed 
can be applied to any quantum systems. Especially, the 
DFF approach can be used to study various quantum 
phase transitions (like the structural phase transition) 
based the LDA calculations. 
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